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ABSTRACT 

The spin-induced quadrupole moment of a rapidly rotating star changes the orbital dy- 
namics in a binary system, giving rise to advance (or regression) of periastron and precession 
of the orbital plane. We show that these effects are important in the recently discovered 
radio pulsar/main sequence star binary system PSR J0045— 7319, and can reliably account 
for the observed peculiar timing residuals. Precise measurements of the apsidal motion and 
orbital plane precession can yield valuable information on the internal structure and rotation 
of the star. The detection of orbital precession implies that the spin of the companion star 
is not aligned with the orbital angular momentum, and suggests that the supernova gave the 
pulsar a kick out of the original orbital plane. Tidal excitation of g-mode oscillations in the 
PSR J0045— 7319 system induces an orbital period change of order |AP or b/P rb| ~ 10~ 6 at 
each periastron passage, but the secular trend depends on the radiative damping time of the 
g-modes. We also discuss the spin-orbit coupling effects for the accreting X-ray pulsars and 
the other known radio pulsar/main sequence binary, PSR B1259— 63. 

Subject headings: binaries: close — pulsars: individual: PSR J0045— 7319, PSR B1259— 63 
— stars: early type — stars: rotation — stars: oscillations 
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1. INTRODUCTION 



The discovery of radio pulsars in binary systems with massive main sequence stars 
(PSR B1259-63, Johnston et al. 1992, 1994 and PSR J0045-7319, Kaspi et al. 1994) makes 
it possible to measure hydrodynamical effects on the pulsar orbit. Indeed, the detection 
of large (±30 msec) frequency-independent timing residuals from PSR J0045— 7319 (Kaspi 
et al. 1995a) has motivated our consideration of the dynamical effects caused by the spin- 
induced quadrupole on the stellar companion. This is expected to be the dominant effect, as 
massive main sequence stars are fast rotators and typically have rotation periods on the order 
of days (Jaschek & Jaschek 1987). We show that the spin-orbit coupling in these systems is 
indeed significant, and leads to both apsidal motion and precession of the orbital plane that 
can account for the timing residuals seen in PSR J0045— 7319. 

The information learned from pulsar / main sequence binaries will complement the studies 
of double main sequence binaries, many of which have large enough quadrupoles (either 
from tides or fast rotation) to cause measurable apsidal motion in their eccentric orbits. As 
reviewed by Claret & Gimenez (1993), the apsidal motion is combined with information on 
the stellar masses, radii and spins to yield the apsidal constant, k, which is a dimensionless 
measure of the density concentration of the stellar interior. These measurements compare 
favorably with those expected from theoretical stellar structure models (see Claret &: Gimenez 
1993 for the few exceptions). The pulsar systems considered here are most sensitive to the 
spin-induced quadrupole, permitting better studies of the orientation and magnitude of the 
spin of the stellar companion. Both of these quantities are important to constrain and/or 
measure in the context of mass transfer prior to the supernova and "kicks" during the neutron 
star birth. 

We start in §2 by explaining why the spin-induced quadrupole is important for these 
systems and then summarize the full equations for the precession and apsidal motion, correct- 
ing some previous typographical errors. We then show in §3 that the induced time delay is 
large and easily measurable for the two radio pulsars under consideration, PSR J0045— 7319 
and PSR B1259— 63 . We demonstrate that this effect is present in the timing data of PSR 
J0045— 7319. In particular, the unusual residuals observed by Kaspi et al. (1995a) are the 
result of fitting a purely Keplerian orbit to irregularly sampled data in which apsidal advance 
and orbital precession are present. Section 4 contains a simple estimate of the dynamical 
tidal effect that might prove important. We show that tidal excitation of high-order gravity 
modes in the companion of PSR J0045— 7319 near periastron can yield significant changes 
in the orbital period and eccentricity which may be measurable. We conclude in §5 by out- 
lining the prospects for constraining the system parameters of PSR J0045— 7319, and the 
possibilities for similar measurements in accreting X-ray pulsars. 



The neutron star's gravitational field couples to both the tide-induced quadrupole and 
the spin-induced quadrupole of the stellar companion. We begin by comparing the relative 
magnitudes of these two effects. The tidal bulge on the stellar companion has a height h ~ 
M p R*/(M c r 3 ), where M c and R c are the companion mass and radius, M p is the pulsar (neu- 
tron star) mass, and r is the instantaneous center of mass separation. The tidal quadrupole is 
then roughly Qtide ~ kM c R^{h/ R c ) ~ kM p R^(R c /r) 3 . We compare this to the quadrupole 
induced by rotation at angular velocity, £l s = 2n/P s , Q sp i n ~ kM c R 2 c (Q? s R 3 / GM C ) . At 
periastron, the ratio of these two is 
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where P OT h and e are the orbital period and eccentricity. Equation (1) shows that when f2 s 
is larger than the orbital angular velocity at periastron, Q or b, P = ^ rb(l + e ) 1 ^ 2 /(l — e ) 3 ^ 2 > 
the spin-induced quadrupole effect is more important than the tidal effect. Binary systems 
containing early-type stars (in particular those with convective cores and radiative envelopes) 
with orbital periods longer than 2 — 3 days are not expected to circularize and synchronize 
appreciably in a time comparable to the main sequence life span, at least when the eccen- 
tricities are not too large (Zahn 1977,1984). Both PSR B1259-63 and PSR J0045-7319 are 
of this type and, as we describe in §3, are timed with sufficient accuracy for these effects to 
be measurable. 

2.1 Spin-Induced Distortion 

The quadrupole from the rotational distortion of the finite-sized companion is specified 
by the difference between the moments of inertia about the spin axis, I3, and an orthogonal 
axis, I\. Claret & Gimenez (1992) calculated the moment of inertia for slowly rotating main 
sequence stars. We incorporate their results by introducing a parameter A, whose value is 
close to unity, so that I 3 = ti.l\M c R 2 (The parameter A also accounts for the increase of 
the equatorial radius of the rotating star, which can be as large as 3R c /2 near the maximum 
rotation rate). The quadrupole from the rotational distortion of the star is proportional to 
the spin squared. Following the conventional definition of the apsidal motion constant k 
(Cowling 1938; Schwarzschild 1958), we write 

h-h = hM c R 2 A (2) 

where tl s is the dimensionless spin Cl s = £} S /(GM C / R^) 1 / 2 of the companion. For stars with 
uniform density, k = 3/4. For the M c ~ 10M Q main-sequence star of interest here, k ~ 0.01, 
depending on the stellar age (Schwarzschild 1958; Claret k, Gimenez 1992). 

2.2 Spin-Induced Apsidal Motion and Orbital Precession 

Because of the spin-induced quadrupole moment, the additional 1/r 3 potential term 
between the components leads to apsidal motion (advance or regression of the longitude of 
periastron) and, when the spin angular momentum S is not aligned with the orbital angular 
momentum L, a precession of the orbital plane. The general expressions for the rates of 
apsidal motion and orbital precession have been derived by Smarr & Blandford (1976) and 
Kopal (1978) using standard celestial mechanics perturbation theory. In the coordinates of 
Figure 1, the rate of change of the dynamical longitude of periastron (measured from the 
ascending node in the invariable plane perpendicular to J = L + S), x, is given by 

where a is the semi-major axis of the (relative) orbit, 9 is the angle between L and S, and 
9 C is the angle between L and J. 2 

2 This can be derived following the perturbation procedure of §11-3. C of Goldstein (1980), 
except that here the action variables are J\ = 27rLcos# c , J2 = 27rL, and note that the 
perturbation Hamiltonian (11-58) depends on 9 (or i in Goldstein's notation), the angle 
between L and S. Note that in Eq. (3.10) of Smarr & Blandford (1976), sin 2 (5 + 9) should 
be sin 2 (6 + 9); also in Eq. (11-60) of Goldstein (1980), the factor (1 - e 2 )" 1 should be 
(1-e 2 )" 2 . 
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The orbital plane precession rate fl prec = <i>J (cf. Fig. 1) can be obtained more directly 
by calculating the interaction torque between the stars. We have 



n _ 3GM P (I 3 - Ji)cosfl J 3tt(J 3 - h) / sin 9 cos 9 \ - 

prec ~ 2a3(l-e 2 )3/ 2 ~SL ~ ~ M c a 2 (l - e 2 ) 2 P mh V sintf, ' ' ; ' 



where the minus sign implies that L precesses in a direction opposite to J, and we have used 
L = [GM 2 M 2 a(l — e 2 )/M t ] 1 / 2 for the Keplerian orbital angular momentum. In the limit of 
small mass ratio M p <C M c and S = |S| S> L = |L|, equations (3) and (4) reduce to the rates 
of apsidal motion and orbital precession for the Earth-satellite system (Goldstein 1980). By 
contrast, for binary pulsars of interest in this paper, L S> S. Using the moment of inertia 
Is = 0.1XM C R 2 the quadrupole from equation (2) we obtain 

10kM n /i?A 3/2 ft, cos 



Q ™ h \(M c M t y/i { a ) (l-e 2 )^/ 2 ' (5) 

where M t = M p + M c . The angle of the precession cone (i.e., the angle between L and J) is 
simply 9 C ~ Ssm6/L in the L 3> S limit. 

The above derivation assumes the companion star behaves as a rigid body. One might be 
concerned that fluid stars respond to external forces differently than a rigid object (Smarr & 
Blandford 1976; Papaloizou & Pringle 1982). This is by no means a resolved issue, although 
it seems unlikely that deviation from rigid body behavior is large since the tidal distortion 
of the star is small compared to the rotational distortion. The observations of these pulsars 
may well provide information on this subtlety. 

2.3 Effects on the Observables 

The apsidal motion and orbital plane precession change u, the observational longitude 
of periastron (measured from the ascending node in the plane of the sky), and i, the orbital 
inclination angle (the angle between L and n, the unit vector along the line-of-sight). Let i 
be the angle between J and n (see Fig. 1), in which case the observational angles u and i 
are related to x an d $ by: 

sin ijj = [(sin 9 C cos i Q + cos 9 C sin i cos <3>) sin x + sin i Q cos x sm ®] > ^ 
cos i = — sin 9 C sin i cos $ + cos i a cos 9 C . 

For L 3> S or 6 C <C 1, these reduce to u> c± x + ^ an d i — io + C cos $. Using equations (2), 
(3) and (4), we then have 

3vr(/ 3 -/i) / 3 . 2 A n kR 2 c Cl 2 s / 3 . 2 A 
U " Af c a 2 (l-e 2 ) 2 P orb i 1 " 2 Sm "J " ^^T^) 2 ( X " 2 Sm ^ ' (?) 



and 

di , 37T(/3-/l) „;„ a „„^„;„^..o 



orb " 



■ sin 6* cos 9 sin $. (8) 



dt M c a 2 (l-e 2 ) 2 P orb U1 a 2 (l-e 2 ) 2 

Thus for typical 9 and <£, the numerical value of di/dt is comparable to that of to. Note 
that the observed apsidal motion is an advance for sin# < (2/3) x / 2 and a regression in the 
opposite case. In addition, twice in one precession cycle, at $ = and it (i.e., when J,L 
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and n lie in the same plane), the rate di/dt for the change of the orbital inclination angle is 
identically zero. 

3. SPIN-ORBIT EFFECTS ON BINARY PULSAR TIMING 

We now consider the effects of spin-orbit coupling on the timing of the two binary pulsar 
systems with main-sequence-star companions. The orbital motion gives rise to a delay of 
T(t) = r p -n/c = r p (t) sin ^(t) sini(t) /c in the pulse arrival time, where r p = rM c /(M p + M c ) 
is the pulsar position vector. The orbital phase (or true anomaly) ^(t) measured from the 
observational ascending node is given by Hf(t) = ^(t)x + &t, where \&(t)x is the Keplerian 
value. The longitude of periastron w(t) and the orbital inclination angle i(t) vary according 
to equations (7) and (8), with the precession phase given by $(t) = O prec t + <J> . 

We are interested in the residual ST(t) = r p -n/c — (r p -n/c)x of the time delay compared 
to the Keplerian value. For t <C l/|fi pre c| and t <C 1/H, we have 6T = 5T aps + 5T prec , where 
the contributions from apsidal motion and orbital plane precession are given by (assuming 
ST = at t = 0) 

r (t) 

ST(t) aps ~ -^sin^cos^i) tbt, 

c (9) 
r (t) di 
ST(t) prec ~ -^i cos i sin (t) — t. 

Clearly, the residual resulting from the neglect of varying uj and i in the timing model 
increases with t, and is modulated by the orbital motion 3 . Equations (9) are the leading 
order corrections for this effect, and are valid for time scales much less than the precession 
period, which is ~ 300 years for PSR J0045— 7319 (see §3.1). The next order terms are of 
order Cjt or (di/dt)t smaller and have different dependences on i and The accumulation 
of accurate timing data over a longer baseline will further constrain the system. 

3.1 PSR J0045-7319 

Radio timing observations by Kaspi et al. (1994) of the 0.93 s pulsar PSR J0045— 7319 
have identified it to be in an eccentric (e = 0.8080) 51.17 day orbit having projected semi- 
major axis OpSini/c = 174 s. They made optical observations in the direction of the pulsar 
and found a B1V star with B = 16.03 and V = 16.19 which showed no evidence for emission 
lines. Combining the effective temperature {T e g ~ 2.4 x 10 4 K) and reddening (Ay ~ 0.27) 
with the observed colors gives R c ~ (6.4 ± 0.5)i? Q (d/59 kpc). The optical radial velocity 
measurements of the B star fix its mass to be (8.8 ± 1.8) Mq for M p = 1.4M© and give 
i = 44 ± 5 degrees (Bell et al. 1995). This then yields a ~ 126^0 ~ 20R C . The rotational 
velocity v rot sm8 sn ~ 113 ± 10 km/s (where 9 sn is the angle between S and n) obtained from 
spectral line broadening (Bell et al. 1995) corresponds to Q s ~ 0.22/ sin^ sn , consistent with 
a rapid rotation near break-up, Cl Sj7nax ~ 0.5. The ratio of f2 s and the angular velocity at 
periastron £l rb, P is ~ 1.1/ sin# sn , confirming that the system is not synchronized. Equations 



3 Note that for very eccentric systems, the changes in to and i mainly occur near periastron, 
rather than accumulating uniformly throughout the orbital period. We neglect this finer 
point in equation (9). 
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(5), (7) and (8) then give 




,0.01 J V a J \ 0.5 / V 2 / y r 

di o / A; \ /20i? c \ Y n s \ . „ „ . ^ rad 
~ 2.3 x 10 -3 — — -4 sin (9 cos 6> sin $ 



(10) 



dt " \0.0l J V a J V 0.5 V yr 

The timing residual AT accumulated in one orbit is estimated by setting r p ~ 2a p and 
* — -P or b/2 in equation (9), giving 

a„sini / fc \{20R c \ 2 fn s \ / 3 . 2 A 

AT prec * °£^P wb ^ 56 (MM Y ^Vf — ainecosesin* ) ms. 
p c V 0.01 A a / I 0.5 i Vsini ' 



(11) 



The modulated residual in just one orbit is thus comparable to the residuals found by Kaspi 
et al. (1995a). 

There are other effects that can also give rise to secular changes in oj and i. The standard 
formula for the apsidal motion due to static tide (Cowling 1938; Schwarzschild 1958) gives 
^tide — 1-3 x 10 _4 (fc/0.01)(20-R c /a) 5 rad/yr, about an order of magnitude smaller than the 
spin-induced Co. This secular trend is degenerate with the spin-induced quadrupole effect. 
Tidal effects do not produce changes in the orbital inclination angle. General relativity also 
induces periastron advance (Landau & Lifshitz 1962) cj gr ~ 6.7 x 10 -5 rad/yr, comparable 
to the tidal effect. The geodetic precession of the companion star's spin (Barker & O'Connell 
1975; note that the spin angular momentum of the pulsar is much smaller than that of the 
companion) results in an orbital precession f2 preC)gr ~ 5.9 x 10~ 6 rad/yr, three orders of 
magnitude smaller than the Newtonian spin- induced quadrupole effect. 

3.2 PSR B1259-63 

This 47.76 ms pulsar is in an eccentric orbit with an Be star companion with P or b = 1236 
days, e = 0.8698, and a p sinz/c = 1296 s (Johnston et al. 1992, 1994). The optical data for 
this emission line star are still inadequate to place strong constraints on the spectral type and 
hence radius of the star, so we will use the fiducial values chosen by Johnston et al. (1994), 
R c = 6R Q and M c = 10M Q . For M p = 1.4M Q , we obtain a = 1091i? o = 182i? c . For 
typical values (A ~ 1, k ~ 0.01 and Cl s ~ 0.5) we get tt prec — —4.1 x 10 -5 cos# rad/yr, 
u ~ 2.4 x 10" 6 (1 - 3/2 sin 2 9) rad/yr, and di/dt ~ 2.4 x 10" 6 sin # cos # sin $ rad/yr. The 
timing residual accumulated in one orbit is given by 

2 ( 12 ) 

AT prec ~ 10 ( I I I I t4 I I -^-sm^cos^sm^ ] ms 



0.0.1 J V a J \ 0.5 / V sin / 
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Thus AT can potentially affect the timing parameters for this system as well. Indeed, using 
a pure Keplerian model for the system resulted in substantial systematic timing residuals 
(Johnston et al. 1994; Manchester et al. 1995). From our estimate above, it is clear that 
these residuals may well have a significant contribution from the unfitted Co and di/dt. A 
full understanding of this system is complicated by the dispersion and radio eclipse caused 
by the wind from the Be star. 

The tide-induced apsidal motion uitide — 2-7 x 10~ 10 (A:/0.01)(182i? c /a) 5 rad/yr, as well 
as the general relativistic effects, are completely negligible for this system. 

Note that the expressions given in §2 describe the secular change of the orbit, i.e., they 
are obtained by orbital averaging. To model closely-sampled timing data spanning only a 
few orbits, as in PSR B1259— 63, a numerical integration of the orbit with the perturbing 
potential may be required. 

3.3 Data Fitting and the Observed Residual in PSR J0045-7319 

Using multi- frequency radio timing observations of PSR J0045— 7319, Kaspi et al. (1995a) 
showed that the observed barycentric pulse arrival times deviated significantly from those 
expected from a pure Keplerian orbit. In particular, they found significant frequency- 
independent timing residuals for a ~10 day duration around periastron. These residuals 
showed systematic variations on time scales of a few days after subtraction of a Keplerian 
orbital fit to the entire 3.3 year data span. The residuals from two periastron observations 
separated by 11 orbits showed very different trends, inconsistent with a simple error in any 
of the Keplerian parameters. They noted that such residuals were unlike those of any other 
pulsar, isolated or binary. 

We show here that the unusual PSR J0045— 7319 residuals can be explained by the 
relatively long-term dynamical effects described above, in spite of the former's apparent 
short-time-scale variations. To do this, we simulated timing observations for a binary pulsar 
which has Keplerian parameters similar to those of PSR J0045— 7319, but in addition has two 
post-Keplerian parameters, the apsidal advance Co and a time-variable projected semi-major 
axis x, where x = a p s'mi. In the simulation, x and uj were incremented by 0.0044 lt-s and 
0.0035° per orbit at periastron, which gave the best match to the observed residuals. These 
numbers correspond to Co ~ 4.4 x 10~ 4 rad/yr and di/dt ~ 1.8 x 10~ 4 rad/yr, consistent 
with equation (16). The actual fit to the data from PSR J0045— 7319 and the implications 
for that source in particular will be discussed elsewhere (Kaspi et al. 1995b). 

The simulated data were then fitted with a pure Keplerian orbit, holding all post- 
Keplerian parameters fixed at zero, thus simulating the analysis done by Kaspi et al. (1995a). 
The residuals for this fit are shown in the upper panel of Figure 2. The effect of the post- 
Keplerian parameters is obvious as an orbital period modulation of increasing amplitude 
relative to the fitting epoch (cf. eq.(9)). The fitted Keplerian parameters are biased away 
from the assumed values because of the least-squares fitting procedure; indeed the fitted 
solution shown in the Figure has "spikes" for this reason. Next, using the same set of 
simulated arrival times, we selected those that coincide with actual observations of PSR 
J0045— 7319 at the Parkes observatory, in order to simulate the sampling interval used by 
Kaspi et al. (1995a) in their timing analysis. Again we fit a simple Keplerian orbit, holding 
all post-Keplerian parameters fixed at zero. The resulting post-fit residuals are shown in the 
middle panel of Figure 2. Close-up views of the two epochs of periastron investigated by 
Kaspi et al. are shown in the bottom panels; a comparison with the data indicates that this 
model closely mimics the observed trends. Thus, we have shown that the unusual residuals 
for PSR J0045— 7319 can be explained as being a result of fitting a purely Keplerian orbit 
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to infrequently and irregularly sampled data in which post-Keplerian apsidal advance and 
orbital plane precession are significant. 



4. POSSIBLE ORBITAL PERIOD CHANGE DUE TO DYNAMICAL TIDE 

The spin-induced quadrupole effects (and the standard tide-induced apsidal motion) 
discussed in the previous sections assume that at a given binary separation the star relaxes 
instantaneously to the equilibrium state. However, dynamical tidal effects associated with 
the excitation of g-modes in the companion star result in energy and angular momentum 
transfer between the star and the orbit during each periastron passage, therefore changing 
in the orbital period and eccentricity. 4 

The energy transfer AE during each periastron passage can be estimated with the 
formalism of Press & Teukolsky (1977), who considered tidal excitations for a parabolic 
orbit. 5 In the quadrupole order, this is given by 



R, 



c 



where r min is the binary separation at periastron, rj = (M c lMt) x ^ 2 {r m i n i 'R c ) 3 ^ 2 is the ratio 
of the time for periastron passage and the stellar dynamical time, and the dimensionless 
function T2{rj) (defined in Press & Teukolsky 1977) involves I = 2 non-radial stellar oscillation 
modes and the orbital trajectory. The resulting fractional change of the orbital size is simply 
Aa/a = —AE/E, and the orbital period change is given by 

|AP orb | _ 3|Aa| M p f R,^ ' 



orb 



2a M r \ a 



3ttM- (l-e)- b T 2 ( V ). (14) 



The angular momentum transferred during periastron is of order AL ~ AE(GM C / R^,)~ 1 / 2 , 
thus \AL/L\ is typically negligible compared to \AE/E\. The resulting change in orbital 
eccentricity is then Ae ~ (Aa/a)(l — e 2 )/(2e). 

To estimate the importance of the dynamical tidal effects, we compare Aa/a with the 
fractional change (due to the static tide) of uj in one orbit Aw/(27r) = u>tide-forb/(27r). For 
PSR J0045— 7319, using the parameters of §3.1, we have rj ~ 7.0, and T 2 (r/) ~ 3.5 x 10~ 4 
for an n = 3 polytrope with adiabatic index T t = 5/3 (e.g., Lee & Ostriker 1986), where 
most contributions to T 2 (ri) come from g-modes of radial order 5 — 9. This gives |Aa/a| ~ 
0.lAw/(27r), implying that the perturbation due to energy transfer is much smaller than the 
tide-induced apsidal motion, and is therefore even smaller than the spin effects discussed in 
§2-3. This estimate agrees with the more detailed study by Kumar, Ao & Quataert (1995). 
The orbital period change due to the energy transfer is given by | AP orD | /P OT h ~ 10 -6 , more 
than an order of magnitude larger than the measurement accuracy of P or b attained over 



4 Mass loss from the companion star can in principle induce similar changes in the orbit. 
However, the absence of eclipses or orbital phase dependent DM variations constrains the 
ionized mass loss rate M c & l(r lo M yr" 1 for PSR J0045-7319 (Kaspi et al. 1994, 1995a), 
suggesting that the dynamical importance of the mass loss is small. 

5 This formalism is also reasonably accurate for highly elliptic orbits. This is because most 
of the energy transfer takes place near periastron, where the kinetic energy and potential 
energy ~ GM p M c /r min are much larger than the binding energy GM p M c / (2a), and the orbit 
resembles a parabolic trajectory. 
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many orbits (Kaspi et al. 1995a). The associated eccentricity change |Ae| ~ 10~ 7 will be 
more difficult to measure. For PSR B1259— 63, again using the parameters of §3.2, we have 
r\ ~ 108, for which T2{rj) is extremely small, and therefore the dynamical tidal effects are 
therefore completely negligible 6 . 

It is important to note that the actual values and signs of AP or b and Ae depend on the 
damping of the tidally excited g-modes. Depending on detailed stellar models, the radiative 
damping time for g-modes in an early-type main sequence star is of order 100 years, although 
some g-modes can be damped in less than a year (e.g., Unno et al. 1979). If the damping time 
£damp is much longer than the orbital period, the energy transfer for an elliptical orbit will 
depend on the phases of the oscillations of different modes, and varies (both in magnitude 
and sign) from one passage to another. In this case, equation (14) only represents the typical 
magnitude of AP or b during each passage, not a steady P OT b', the averaged long-term orbital 
decay rate is expected to be of order — |AP or b|/tdamp- On the other hand, if the damping time 
tdamp is much shorter than the orbital period, the energy transfer near periastron is a one-way 
process, i.e., only from the orbit to the star. Therefore AP or b/-P rb and Ae represent steady 
negative values, i.e., the orbit decays and circularizes. If such a short damping timescale 
indeed applies to the high order g-modes excited in PSR J0045— 7319, the orbital decay can 
be observable in the timing data. A more accurate calculation of this effect is complicated 
by the rapid rotation of the B-star, which changes the g-mode structure when the rotation 
frequency is comparable to the g-mode frequencies. 

5. DISCUSSION 

Smarr & Blandford (1976) first discussed the Newtonian spin-orbit effects in the con- 
text of PSR B1913+16. That companion turned out to be a neutron star, which renders 
these Newtonian effects completely negligible and facilitates excellent tests of general relativ- 
ity (e.g., Taylor & Weisberg 1989). The recent discovery of two radio pulsar/main sequence 
star binaries, PSR J0045— 7319 and PSR B1259— 63, finally allows for measurements of hy- 
drodynamical effects in binary pulsar systems. In this paper we have shown that the Newto- 
nian spin-orbit coupling is particularly important in PSR J0045— 7319 and PSR B1259— 63 
systems. 

5.1 Constraints on the Binary Systems and Evolution 

When the spin-induced quadrupole is the dominant effect, measuring Co and di/dt pro- 
vides valuable information on the system: (i) the sign of Co constrains the spin-orbit angle 9 
(cf. eq.(7)); (ii) the ratio of Co and di/dt yields a relation between 9 and $; (iii) the angle 9 sn 
is a function of 9, and i. Knowing i and by choosing a reasonable range of 9 sn (e.g., for the 
PSR J0045— 7319 system: since v ro t must be smaller than the break-up limit ~ 512 km/s, we 
have | sin# sn | ^ 113/512 = 0.22), we obtain another relation between 9 and (iv) With the 
angles determined from (i)-(iii), the observed values of Co and di/dt then constrain the stellar 
structure and rotation rate of the companion. In addition, long-term timing observations 
will measure various angles in addition to Co and di/dt, yielding additional constraints on the 
system parameters (see discussion following eq. (9)). 

Other effects may prove important for these systems, including the apsidal motion due 
to static tide and general relativity. The dynamical tide (§4) may induce measurable changes 

6 Kochanek (1993) has previously discussed the dynamical tidal effects in the context of 
PSR B1259— 63. However, the observational data available then did not allow him to draw 
any firm conclusion. 
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in P or b and e. However, a detailed theoretical treatment of this effect is complicated by the 
rapid rotation of the companion star. 

The information obtained from the systems can have important implications for the 
evolution of neutron star binaries. For example, if mass transfer in the pre-supernova binary 
forces S to be parallel to L (a very likely scenario), then a non-zero spin-orbit angle in 
the pulsar/main-sequence star system implies that the supernova gave the neutron star a 
kick velocity out of the original orbital plane. From the values of 9 and e of the current 
system (assuming that circularization and spin-orbit alignment are not efficient in the current 
system), one may actually constrain the kick velocity. Cordes & Wasserman (1984) discussed 
how such kicks might misalign the neutron star spin from the orbital angular momentum in 
PSR 1913+16, giving rise to measurable geodetic precession. 

5.2 X-ray Binaries 

Insight into orbital plane precession might prove important for accreting X-ray pulsars as 
well. Apsidal motion has been searched for in eccentric accreting systems. In particular, the 
2a upper limit on the apsidal motion in the eccentric 8.9 day binary Vela X-l (which contains 
a supergiant) is Co = 1.6 deg/yr (Deeter et al. 1987). Tamura et al. (1992) claimed a detection 
of apsidal motion in the eccentric Be system 4U 0115+63 at the level Co = 0.030+0.016 deg/yr. 
More recent observations of this source with the BATSE instrument (Cominsky, Roberts &: 
Finger 1994) contradict the claim of Tamura et al. (1992), find no evidence for apsidal motion. 
However, these measurements could be biased by the fact that a x sin? was assumed constant. 

The less accurate arrival times makes the search for orbital plane precession in the ac- 
creting X-ray pulsars more difficult than for radio pulsars. However, the orbital parameters 
are sufficiently accurate that one can search for orbital precession through temporally sep- 
arated measurements of a x sini, the projected size of the orbit for the X-ray pulsar. The 
system 4U 0115+63 with P or b = 24 days and e = 0.34 may experience Ai ~ 6 x 10~ 3 rad 
over a 17 year baseline (from equation (8) assuming a 1OM main sequence star with radius 
R = 6Rq, k = 0.01 and rotating slowly enough, fl s ~ 0.05, or accidentally having sin 2 9 
close to 2/3 so as not to violate the upper limit on Co from Cominsky et al. 1994). This is 
a factor of 6 larger than the fractional error on a x s'mi measured in 1978 (Rappaport et al. 
1978) and may well allow for a very sensitive measurement of the rotation of the Be star in 
this system. 

We thank Roger Blandford and Sterl Phinney for helpful discussion. This work has 
been supported in part by NASA Grants NAGW-2394 to Caltech and NAG5-2819 to U.C. 
Berkeley. L. B. was supported as a Compton Fellow at Caltech through NASA grant NAG5- 
2666 during part of this work. V.M.K. thanks Dave Van Buren and is supported by a 
Hubble Fellowship through grant number HF-1061.01-04A from the Space Telescope Science 
Institute, which is operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS5-26555. Research by V.M.K. was carried out at the Jet 
Propulsion Laboratory, California Institute of Technology, under a contract with the National 
Aeronautics and Space Administration. 



10 



REFERENCES 



Barker, B. M., & O'Connell, R. F. 1975, ApJ, 199, L25 

Bell, J. F., Bessell, M. S., Stappers, B. W., Bailes, M. & Kaspi, V. 1995, ApJL, in press 
Claret, A., & Gimenez, A. 1992, A&A Suppl., 96, 255 
. 1993, A&A, 277, 487 

Cominsky, L., Roberts, M. & Finger, M. H. 1994 in The Second Compton Symposium, ed. 

C. E. Fichtel, N. Gehrels, & J. P. Norris (AIP: New York) p. 294 
Cordes, J. M. & Wasserman, I. 1984, ApJ, 279, 798 
Cowling, T. G. 1938, MNRAS, 98, 734 

Deeter, J. E., Boynton, P. E., Lamb, F. K., & Zylstra, G. 1987, ApJ, 314, 634 
Goldstein, H. 1980, Classical Mechanics, 2nd Ed. (Addison- Wesley), p. 512 
Jaschek, C, & Jaschek, M. 1987, The Classification of Stars (Cambridge Univ. Press: New 
York) 

Johnston, S. et al. 1992, ApJ, 387, L37 
Johnston, S. et al. 1994, MNRAS, 268, 430 
Kaspi, V. M., et al. 1994, ApJ, 423, L43 

Kaspi, V. M., Manchester, R. N., Bailes, M. & Bell, J. F. 1995a, Compact Stars in Binaries, 
IAU Symp. 165, ed. J. van Paradijs, E. P. J. van den Heuvel, & E. Kuulkers (Kluwer 
Dordrecht), in press. 

Kaspi, V. M., et al. 1995b, in preparation. 

Kochanek, C. S. 1993, ApJ, 406, 638 

Kopal, Z. 1978, Dynamics of Close Binary Systems (D. Reidel Pub.: Dordreicht, Holland), 
p 232 

Kumar, P., Ao, C. O., & Quataert, E. J. 1995, ApJ, in press 

Landau, L. D., & Lifshitz, E. M. 1962, Classical Theory of Fields (Oxford: Pergamon) 

Lee, H. M., & Ostriker, J. P. 1986, ApJ, 310, 176 

Manchester, R. N., et al. 1995, ApJL, in press 

Papaloizou, J. C. B., & Pringle, J. E. 1982, MNRAS, 200, 49 

Press, W. FL, & Teukolsky, S. A. 1977, ApJ, 213, 183 

Rappaport, S., Clark, G. W., Cominsky, L., Joss, P. C. & Li, F. 1978, ApJ, 224, LI 
Schwarzschild, M. 1958, Structure and Evolution of the Stars (Princeton University Press) 
Smarr, L. L., & Blandford, R. D. 1976, ApJ, 207, 574 

Tamura, K., Tsunemi, H., Kitamoto, S., Hayashida, K. & Nagase, F. 1992, ApJ, 389, 676 
Taylor, J. H. & Weisberg, J. M. 1989, ApJ, 345, 434 

Unno, W., Osaki, Y., Ando, FL, & Shibahashi, H. 1979, Nonradial Oscillations of Stars 

(Univ. of Tokyo Press), pp. 205-217. 
Zahn, J.-P. 1977, A&A, 57, 383 

. 1984, in Observational Tests of the Stellar Evolution Theory, eds. A. Maeder 

and A. Renzini (Dordrecht: Reidel), p. 379. 



11 



Figure Captions 

FIG. 1. — Binary geometry and the definitions of different angles. The invariable plane 
(XY) is perpendicular to the total angular momentum vector J = L + S, and the line-of- 
sight unit vector n lies in the YZ plane, making an angle i a with J. The inclination of the 
orbit with respect to the invariable plane is 6 C (which is also the precession cone angle of L 
around J), while the orbital inclination with respect to the plane of the sky is i. The angle 
between L and S is 6. The orbital plane intersects the invariable plane at ascending node 
A, with a longitude $ measured in the invariable plane (<I> is also the phase of the orbital 
plane precession). The dynamical longitude of periastron (point P), measured in the orbital 
plane, is x- The observational longitude of periastron is lo (not shown in the Figure). 

FIG. 2. — Simulated pulse arrival times for PSR J0045— 7319. The upper panel shows the 
residuals from a fit of a pure Keplerian orbit to fake pulse arrival time data ("measured" 
daily for four years) for a binary pulsar like PSR J0045— 7319 in which apsidal advance and 
orbital plane precession are important. The middle panel shows the residuals from a fit 
of this fake data sampled only at those epochs coinciding with actual observations of the 
pulsar. The two vertical solid lines show epochs of periastron separated by 11 orbits. The 
lower panels show close-ups of those epochs; the observed trends in the residuals are similar 
to those observed by Kaspi et al. (1995a). 
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